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Abstract. 

By a modification of the method that was applied in (Korolev and Shevtsova, 2009), here 
the inequalities 

,„ ^ 0.335789(/? 3 + 0.425) 

and 

,„ ^ 0.3051(/? 3 + 1) 



n 



are proved for the uniform distance p(F n ,$>) between the standard normal distribution function 
<J? and the distribution function F n of the normalized sum of an arbitrary number n ^ 1 of 
independent identically distributed random variables with zero mean, unit variance and finite 
third absolute moment /3 3 . The first of these inequalities sharpens the best known version of the 
classical Berry-Esseen inequality since 0.335789(/3 3 + 0.425) < 0.335789(1 + 0.425)/? 3 < 0.4785/3 3 
by virtue of the condition (3 3 ^ 1, and 0.4785 is the best known upper estimate of the absolute 
constant in the classical Berry-Esseen inequality The second inequality is applied to lowering 
the upper estimate of the absolute constant in the analog of the Berry-Esseen inequality for 
Poisson random sums to 0.3051 which is strictly less than the least possible value of the absolute 
constant in the classical Berry-Esseen inequality. As a corollary, the estimates of the rate of 
convergence in limit theorems for compound mixed Poisson distributions are refined. 

Key words: Central limit theorem, Berry-Esseen inequality, smoothing inequality, Poisson 
random sum, mixed Poisson distribution 



*Research supported by the Russian Foundation for Basic Research, projects 08-01-00563, 08-01-00567, 
08-07-00152 and 09-07-12032-ofi-m, and also by the Agency for Education of Russian Federation, state 
contracts P-1181 and P-958. 

^Department of Mathematical Statistics, Faculty of Computational Mathematics and Cybernetics, 
Moscow State University; Institute for Informatics Problems, Russian Academy of Sciences, 
vkorole v@ comt v . ru 

■^Department of Mathematical Statistics, Faculty of Computational Mathematics and Cybernetics, 
Moscow State University, ishevtsova@cs.msu.su 



1 



1 Introduction 



By JF 3 we will denote the set of distribution functions with zero first moment, unit second 
moment and finite third absolute moment (3 3 . Let X 1; X 2 ,... be independent random 
variables with common distribution function F G JF 3 defined on a probability space 
(n,A, P). Denote 



= F* n (xV^) = P ( Xl + --_ +Xn < x 



= / <j,(t)dt, 4>(x) = -^ e ~ x2/2 , x G 



'2tt 

The classical Berry-Esseen theorem states that there exists a finite positive absolute 
constant C which guarantees the validity of the inequality 

p(F n ,$) = sup|F n (x)-$(x)KCo4^ (1) 

In 



for all n ^ 1 and any F G JF 3 (Berry, 1941), (Esseen, 1942). The problem of establishing 
the best value of the constant Co in inequality (1) is very important from the point of view 
of practical estimation of the accuracy of the normal approximation for the distribution 
functions of random variables which may be assumed to have the structure of a sum of 
independent random summands. 

This problem has a long history and is very rich in deep and interesting results. Upper 
estimates for Co were considered in very many papers. Here we will not repeat a detailed 
history of the efforts to lower the upper estimates of C from the original works of A. 
Berry (Berry, 1941) and C.-G. Esseen (Esseen, 1942) to the papers of I. S. Shiganov 
(Shiganov, 1982), (Shiganov, 1986) presented in (Korolev and Shevtsova, 2009). We will 
restrict ourselves only to an outline of the recent history of the subject. 

After some lull that lasted more than twenty years, recently the interest to the problem 
of improving the Berry-Esseen inequality rose again and resulted in very interesting and 
in some sense path-clearing works. In 2006 I. G. Shevtsova improved Shiganov's upper 
estimate by approximately 0.06 and obtained the estimate C ^ 0.7056 (Shevtsova, 
2006). In 2008 she sharpened this estimate to C < 0.7005 (Shevtsova, 2008). In 2009 
the competition for improving the constant became especially keen. On 8 June, 2009 I. 
S. Tyurin submitted his paper (Tyurin, 2009a) to the «Theory of Probability and Its 
Applications». That paper, along with other results, contained the estimate Co ^ 0.5894. 
Two days later the summary of those results was submitted to «Doklady Akademii Nauk» 
(translated into English as «Doklady Mathematics») (Tyurin, 2009b). Independently, on 
14 September, 2009 V. Yu. Korolev and I. G. Shevtsova submitted their paper (Korolev 
and Shevtsova, 2009) to the «Theory of Probability and Its Applications*. In that paper 
the inequality 

,„ ^ x 0.34445(/3 3 + 0.489) 

p{F n ,$)< n^l, (2) 

in 



was proved which holds for any distribution F G T<$ yielding the estimate C ^ 0.5129 by 
virtue of the condition j3 3 ^ 1. Finally, on 17 November, 2009 the paper (Tyurin, 2009c) 
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was submitted to the «Russian Mathematical Surveys» (its English version (Tyurin, 
2009d) appeared on 3 December, 2009 on arXiv:0912.0726vl). In this paper the estimate 
Co ^ 0.4785 is proved. So, the best known upper estimate of the absolute constant Co in 
the classical Berry-Esseen inequality (1) is C ^ 0.4785 (Tyurin, 2009c). 
On the other hand, in 1956 C.-G. Esseen showed that C ^ Ce where 



10 + 3 

C E = v Z. = 0.409732... 

6^27T 



(Esseen, 1956). In 1967 V. M. Zolotarev put forward the hypothesis that in (1) Co = Ce 
(Zolotarev, 1967a), (Zolotarev, 1967b). However, up till now this hypothesis has been 
neither proved nor rejected. 

To prove (2) we used an observation that from inequality (1) it obviously follows that 
for any k ^ there exists a finite positive absolute constant Ck which guarantees the 
validity of the inequality 

p^lXC^ (3) 
Jn 



for all n ^ 1 and F e JF 3 (for example, inequality (3) trivially holds with Ck = C ). 

Following the lines of the reasoning we used in (Korolev and Shevtsova, 2009) to 
prove (2), with the only change in the way of estimation of the difference between 
characteristic functions in the neighborhood of zero (see lemma 2 below), in this paper 
we will demonstrate a special method of numerical estimation of Ck in (3). This method 
yields two special values of k: k = k and k — 1. The first value, /c , minimizes the upper 
estimate of Ck(l + k) yielding the best (within the method under consideration) upper 
estimate of Co in (1) since 

C ^ minC fc (l + k) 

k^0 

by virue of the condition [3 3 ^ 1. At the same time the second value, k = 1, minimizes 
Ck in (3). As we will see, k = 1 plays the main role in improving the absolute constant in 
the analog of the Berry-Esseen inequality for Poisson and mixed Poisson random sums. 
Inequality (3) with k = k and k = 1 is an improvement of the inequality 

B 3 + 1 
p(F n , $) ^ 0.3450- 



n 

we proved in (Korolev and Shevtsova, 2010a). In (Korolev and Shevtsova, 2010b) this 
inequality was applied to sharpening the analog of the Berry-Esseen inequality for Poisson 
random sums and it was for the first time demonstrated that the absolute constant in 
this analog can be made strictly less than that in the classical Berry-Esseen inequality. 

In the papers (Shevtsova, 2010a) and (Korolev and Shevtsova, 2010a) it was shown 
that the constant Ck in (3) cannot be made less than the so-called lower asymptotically 
exact constant in the central limit theorem, that is, 



C k ^ — = 0.2659... 
3V2T 



so that the gaps between the least possible value of the constant Ck and its upper 
estimates given in theorems 1 and 2 below are rather small and do not exceed 0.07 and 
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0.035, respectively, which is important from the point of view of practical applications of 
inequalities (6) and (7). 

Our investigations were to a great extent motivated by a series of results of Hakan 
Prawitz and Vladimir Zolotarev outlined below. 

First, since estimates of the accuracy of the normal approximation for distributions of 
sums of independent random variables are traditionally constructed with the use of the 
so-called smoothing inequalities which estimate the (uniform) distance between the pre- 
limit distribution function of the standardized sum of independent random variables and 
the limit standard normal distribution function through some integral of the (weighted) 
absolute value of the difference between the corresponding characteristic functions, the 
shape of the dependence of the final estimate on the moments of summands is fully 
determined by the shape of dependence of the majorant of characteristic functions on 
these moments. In (Prawitz, 1973) the following result was presented. Let f(t) be the 
characteristic function corresponding to the distribution function F G JF 3 . Denote 

I cos x — 1 + x 2 12\ 
H = sup - = — = 0.09916191... 

x>0 X 6 

and let 9 = 3.99589567... be the unique root of the equation 

3(1 - cos 9) -9 sin 9 - 9 2 /2 = 0, 
lying in the interval (tt, 2tt). Then 



+ + 1)|*|- l*K (/33 + 1) , 

1 - cos ((#, + l)t) 

1 Tfa+W ' *o< (ft + 1)1*1 <2tt, 

2tt 

1, t > 



(A 



As is easily seen, the majorant for \f(t)\ established by this inequality depends on /3 3 
through the function ijj{f3z) — 03 + 1- This is the first hint at that the final estimate for 
p(F n , <3>) should also depend on f3 3 through the function ^(^3) = ^3 + 1. 

Second, in (Prawitz, 1975b) H. Prawitz announced an inequality with unusual 
structure 

tr , ^ 2 03 1 c l (8 3 ) 2 + c 2 03 + c 3 „ ^ 

3V27T y/n-1 2-^2^(71-1) n-1 

where c\, C2 and C3 are some finite positive constants. In the same paper he suggested 
that the coefficient 

= 0.2659... 



3V27T 

at Pz/\/n — 1 cannot be made smaller. Probably, H. Prawitz intended to publish the strict 
proof of (4) in the second part of his work which, unfortunately, for some reasons remained 
unpublished (the title of (Prawitz, 1975b) contains the Roman number I indicating the 
assumed continuation). 
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This Prawitz' inequality (4) seemed to have bepuzzled some specialists in limit 
theorems of probability theory. In particular, it was bypassed in the well-known books 
(Petrov, 1987), (Zolotarev, 1997) (in both of these books there is even no reference to 
any of Prawitz' works). Only in the book (Petrov, 1995) there appears a reference to 
the paper (Prawitz, 1975a) dealing with some estimates for characteristic functions, but 
the paper (Prawitz, 1975b) containing inequality (4) is again ignored. In Mathematical 
Reviews (Dunnage, 1977) there is only a fuzzy remark concerning «some improvements 
for identically distributed summands». Probably, this attitude of some specialists to 
inequality (4) is caused by that at first sight this inequality contradicts the Esseen's 
result that Cq ^ Ce cited above, since 



10 + 3 

< 



3V2vr 6V27T 

However, a thorough analysis of the published part of Prawitz' work convinces that 
inequality (4) is valid. A strict proof of a similar inequality for not necessarily identically 
distributed summands with the third term being 0((/3 3 / v^) 5//3 ) was given by V. Bentkus 
(Bentkus, 1991), (Bentkus, 1994) (for identically distributed summands, the result of 
Benkus is slightly worse than (4) where the third term is O '((/^ / V™) 2 )) • 

Inequality (4) has a very interesting structure: from the main term of order 0(n -1 / 2 ) 
of the estimate of the accuracy of the normal approximation a summand of the form 
l/y/n is separated. This summand may be considerably less than the Lyapunov fraction 
03/Vn. Moreover, in the double array scheme it may happen so that even if the Lyapunov 
condition (3^/y/n — > holds, the quantity (3% = fain) may infinitely increase as n — > oo 
so that the summand of the form n~ 1//2 is infinitesimal with a higher order of smallness 
than the Lyapunov fraction (3 3 (n)/^/n. Thus, inequality (4) is the second hint at that in 
a reasonable estimate of p(F n , <3>) depending on /5 3 the term of order 0(n~ 1 / 2 ) should be 
split into two summands of the form fe/y/n and 1/y/n respectively. 

By the way, speaking of the history of inequality (4), it has to be noted that actually 
it is a further development of the inequality 

,„ _ 0.32/53 + 0.25 

p(F n ,$K , n^3, (5) 

\/n — 2 



which holds under the condition \/n — 1 ^ 3.9(^3 + 1). The proof of (5) was given by 
H. Prawitz in his lecture on 16 June, 1972 at the Summer School of the Swedish Statistical 
Society in Lottorp (Prawitz, 1972a). 

So, the final shape of inequality (3) was prompted by the works of H. Prawitz 
mentioned above. As this is so, the main role goes to the problem of a proper estimation 
of the constant C^. To solve this problem we use a method which is a further development 
of the ideas of V. Zolotarev presented in (Zolotarev, 1965), (Zolotarev, 1966), (Zolotarev, 
1967a) and (Zolotarev, 1967b). This method will be described in detail below. 

The paper is organized as follows. In Section 2 the basic results are proved. Namely, 
here we prove inequality (3) with k = k = 0.425 (theorem 1) and with k — 1 (theorem 2). 
In Section 3 theorem 2 is applied to sharpening the analog of the Berry-Esseen inequality 
for Poisson random sums. We show that despite a prevalent opinion that the absolute 



5 



constant in this inequality should not be less than the absolute constant in the classical 
Berry-Esseen inequality, as a matter of fact this is not so and the constant in the Berry 
Esseen inequality for Poisson random sums does not exceed 0.3051, which is, as it has been 
already mentioned, strictly less than the least possible value Ce ~ 0.4097 of the constant 
C in (1). Finally, in Sections 4 and 5 the result of Section 3 is used for improving the 
estimates of the rate of convergence of compound mixed Poisson distributions with zero 
and non-zero means to scale and location mixtures of normal laws, respectively. 

2 The basic results 

2.1 Formulations and discussion 

Practical calculations show that under the algorithm we use for the estimation of Ck (see 
Section 3) the resulting majorant of the constant Ck decreases as k increases from to 1. 
At the same time for ^ k ^ k « 0.425 the obtained estimates of CV(1 + k) remain 
constant, and for k > k they begin to increase although in the interval k < k < 1 the 
obtained estimate of Ck decreases. Thus, we can present two computationally optimal 
values of k in (3): ko = 0.425 and k\ = 1. The first of them delivers the minimum value 
to the upper estimate of Ck(l + k), thus solving the problem of estimation of Co in (1), 
whereas the second, maximin, minimizes the estimate of Ck in (3). 
The use of k = k in (3) gives the following result. 

Theorem 1. For all n ^ 1 and all distributions with zero mean, unit variance and 
finite third absolute moment (3 3 we have the inequality 

,„ _ 0.335789(/3 3 + 0.425) 
p(F n , $ ^ ^= (6) 

Remark 1. Under the conditions imposed on the moments of the random variable X% 
we always have (3 3 ^ 1. Therefore, 

0.335789(/3 3 + 0.425) ^ 0.335789(1 + 0.425)/? 3 < 0.4785/? 3 . 

Hence, inequality (6) is always sharper than the classical Berry-Esseen inequality (1) with 
the best known constant Co = 0.4785 for all possible values of j3 3 , although the same prior 
information concerning the distribution F is required for its validity (namely, only the 
value of the third absolute moment /3 3 ). 

REMARK 2. Inequality (6) is an «unconditional» variant of the «conditional» Prawitz 
inequality (5) and is a practically computable analog of inequality (4) with a slightly 
(approximately by 0.07) worse first coefficient and a slightly better (approximately by 
0.05) second coefficient, but without the third summand that contains unknown constants. 

Remark 3. Even if the hypothesis of V. M. Zolotarev that C = C E = 0.4097... in 
(1) (see (Zolotarev, 1967a), (Zolotarev, 1967b)) turns out to be true, then, due to that 
(3 3 ^ 1, inequality (6) will be sharper than the classical Berry-Esseen inequality (1) for 
f3 3 ^ 1.93. 

The use of k = 1 in (3) yields the following result. 
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Theorem 2. For all n ^ 1 and a/i distributions with zero mean, unit variance and 
finite third absolute moment f3 3 we have the inequality 

REMARK 4. Inequality (7) is another unconditional » variant of the «conditional» 
Prawitz inequality (5). Moreover, the first coefficient in (7) is less than that in (5) by 
approximately 0.02 whereas the second coefficient in (7) is greater than that in (5) by 
approximately 0.05. 



2.2 Proofs of basic results 
2.2.1 Auxiliary statements 

As we have already mentioned above, to prove theorem 1 we will follow the lines of the 
approach proposed and developed by V. M. Zolotarev in his works (Zolotarev, 1965), 
(Zolotarev, 1966) and (Zolotarev, 1967). This approach is based on the application 
of smoothing inequalities which make it possible to estimate the distance between 
distribution functions via the distances between the corresponding characteristic 
functions. Within this approach the key points are: (i) the choice of a proper smoothing 
inequality; (ii) the choice of a proper smoothing kernel in a smoothing inequality; (iii) 
the choice of proper estimates for the distance between characteristic functions; (iv) the 
choice of a proper computational optimization procedure. 

We will describe these points one after another as they are used in the proof of 
theorems 1 and 2. The corresponding statements will have the form of lemmas. 

We begin with the smoothing inequality. In most papers dealing with the estimation 
of the constant in the Berry-Esseen inequality (1) smoothing inequalities of the same 
type were used. This type of smoothing inequalities was introduced by V. M. Zolotarev. 
In the original paper (Zolotarev, 1965), just as in similar inequalities in the earlier 
papers of Berry (Berry, 1941) and Esseen (Esseen, 1942), the kernel was used which 
had a probabilistic sense, that is, which was the probability density of some symmetric 
probability distribution. In the paper of Van Beek (Van Beek, 1972) it was noticed that 
this condition is not crucial. Van Beek proposed to use symmetric kernels with alternating 
signs. Concurrently with (Van Beek, 1972), the paper of V. Paulauskas (Paulauskas, 1971) 
was published in which the original smoothing inequality of Zolotarev was generalized 
(and hence, sharpened) to the case of positive non-symmetric kernels. It is interesting to 
notice that although in the final part of the paper of Paulauskas it was noted that the 
smoothing inequality proved in that paper was destined, in the first place, for improving 
the constant in the Berry-Esseen inequality, as far as we know, unfortunately no one ever 
used the Paulauskas inequality for this purpose. In (Shevtsova, 2009b) a new smoothing 
inequality was proved which generalizes (and hence, sharpens) both Paulauskas' and Van 
Beek's inequalities to the case of non-symmetric kernels with alternating signs. However, 
all these inequalities yield worse estimates than the Prawitz smoothing inequality proved 
in (Prawitz, 1972b). 
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The characteristic function of the standardized sum (Xi+. . .+X n )/ \fn will be denoted 

fn(t), 

oo 

fn(t) = [ e Ux dF n (x), t e E. 



Also denote 



r n (t) = \Ut)-e~ t2 l 2 \. 



LEMMA 1 (Prawitz, 1972b). For an arbitrary distribution function F and n ^ 1 for 
any < to ^ 1 and T > we have the inequality 

p(F n ,$) ^ 2 [° \K(t)\r n (Tt)dt + 2 [ \K(t)\ ■ \f n (Tt)\dt+ 

JO Jtn 



to 



to 



7l 



where 



m = \(i-\t\) + l 



[1 - |t|)cot7rt + 



to 
signi 

71 



-1 < t < 1. 



Remark 5. In (Vaaler, 1985) a proof of a result similar to the Prawitz inequality 
stated by lemma 1 was given by a techniques different from that used in (Prawitz, 1972b) 
and it was also proved that the kernel K(t) defined by (8) is in some sense optimal. 

Now consider the estimates of the characteristic functions appearing in lemma 1. For 

e > set 

V/2-xe|t| 3 , \t\^6 /e, 

X (*,e) = < 1 ~™ Et , e <e\t\^2n, (9) 
e 2 

0, \t\ > 2tx/s, 

where 9 = 3.99589567... is the unique root of the equation 

6 2 + 2#sin# + 6(cosfl - 1) = 0, tt ^ 6 < 2tt, 



x = sup 

x>0 



cosx — 1 + x 2 /2| cosx— l + x 2 /2 



x=e 



(10) 

0.09916191... (11) 



It can easily be made sure that the function x(t, e) monotonically decreases in e > for 
any fixed tel. 

The Lyapunov fraction will be denoted £ = (3 3 / ^/n. In addition, denote 



Lemma 2. For any F £ JF 3; n ^ 1 and t £ E £/ie following estimates take place: 



|/„(t)|< l-- X (t,Q 



n/2 



/i(t,4,n) 



\fn(t)\ ^ex V {- X (tJ n )} = f 2 (t,Q, 



\fn(t)\ ^exp 



Remark 6. Apparently, the function fi(t,£ n ,n) was used in the problem of numerical 
evaluation of the absolute constants in the estimates of the accuracy of the normal 
approximation for the first time in (Korolev and Shevtsova, 2009). The second and the 
third estimates presented in lemma 2 are due to H. Prawitz (Prawitz, 1973), (Prawitz, 
1975b). 

Remark 7. Evidently, fi(t,e,n) ^ f2(t,e) for all n ^ 1, e > and t 6 R. Moreover, 
from the result of Prawitz (Prawitz, 1973) it follows that f2{t,e) ^ fs(t,£) for all e > 
and t e R, thus the sharpest estimate for \ f n (t)\ is given by f\(t, £ n , n), while the estimates 
fj(t, £ n ), j = 2, 3, possess a useful property of monotonicity in £ n which is very important 
for the computational procedure. 

Lemma 3 (Tyurin, 2009a), (Tyurin, 2009c), (Tyurin, 2009d). For any Fef 3 ,01 
and t G R we have 



r n {t) < 



-t 2 /2 



1*1 u 2 



3 u 2 /2 



/ 



U 



n-1 



The combination of lemmas 2 and 3 allows to obtain an estimate for the difference 
of the characteristic functions in the neighborhood of zero, which is sharper than all the 
analogous estimates used in the preceding works: 



r n (t) < £e~^ [ lt %e^ 
Jo * 



2 / „ 1 

i--x(u,e+-= 

n V Jn 



(n-l)/2 



du = n(t,£,n), tel. 



From what was said above it follows that the substitution of the functions fj(t,£ n ), j = 
2, 3, instead of f\(t, £ n , n) into the right-hand side of the last inequality does not make the 
resulting estimate less, thus, we obtain two more estimates for r n {t) which monotonically 
increase in £: 



r n (t)^£e- t ' 2 / 2 /Ve^exp 
o 1 

1*1 „2 



n — 1 / 

X ( u, £ + —= ) \du = r 2 (t, £, n) 



/rl 2 2 

r n (t) ^ £e~ t2/2 J y exp jx£ n w 3 + ^-(l - 2x£ n uj^du = r 3 (t,£,n), te 



(recall that £ n = £ + l/y/n). 
Noticing that 



K(t)- 



2nt 



1 \2 

-), 0<t<l, 



^(1 — y 1 + (cotvrf 

we can estimate p(F n , $) for any n ^ 2 and F with a fixed Lyapunov fraction £ as 
p(F n ,$)^2 P \K(t)\-r 1 (Tt,£,n)dt + 2 [ \K{t)\ ■ fr(Tt, £ + 1/y/E, n) dt+ 

J0 Jtn 
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+1 J°° e -T>e/2dt + £ {1 _ t) ^ 1 + (^ cot7rt _l_y e ~T^/2 dt = D ^ nMiT) 

with arbitrary positive T and to ^ 1- 

The following lemma makes it possible to bound above the set of the values of n under 
consideration when estimating the constant Ck in inequality (3) with < k ^ 1. 

Lemma 4. For any positive N, k ^ 1 anc? £ > k/y/N for all t e K t/ie following 
estimates hold: 



sup (i, 5 + ^-^) = (i, e + ^7=^) = /j,jv(^ e), j = 1, 2 



l-fc \ _ / 1- At 

fc \ _, 2/2 r 1 * 1 u 2 r u 2 N - 1 / 1-A;\ . 
supr 2 (t,e ^,n)<£e 7 / — exp \ — —x[u,e^ j=- ) \du = r 2 ,jv(t, e), 

For 



Jo 2 r l2 JV ' 
t\ < T(iV,£) =min | A^ 1 /V 1/2 ,(2x£)" 1 } 



w;e a/so /iawe i/ie estimate 



sup rs (f , e - -L < ^ f e "' <|3 - e ~* 2/2 = r s (*, e). 

Proof. The first two statements are trivial consequences of the monotonicity of the 
functions x(t,£+ (1 — k)/y/n) and fj(t,e+ (1 — fc)/\/n), j = 1,2, with respect to n ^ 1. 
To prove the third statement note that the function r 3 can be written in the form 

1*1 „2 



where 



r 3 ^,£ — — |=,nj = e * 2 / 2 y ^- exp {xeu 3 + g(n, m)}c?m, 

/ \ n / 1 \ a(u) . . u 2 f \ 

g(x, u) = In ye — — H , x > 0, a(u) = — yl — 2xeuj , u > 0. 

Since |t| ^ {2xe)~ l under the conditions of the lemma, we have a(u) ^ for all u ^ \t\. 
Let us establish that g(x,u) monotonically increases in x ^ N and u ^ T(N,e). Indeed, 
the derivative 

dg(x,u) 1 a(u) 

dx 2x(s\/x — 1) x 2 

is non-negative if and only if x — 2a(u)e\/x + 2a(u) ^ 0. Since a(u) ^ 0, the last condition 
is satisfied, if \fx ^ 2a{u)e = eu 2 {\ — 2xeu), or, particulary, if \fx ^ eu 2 . So, for &\\x ^ N 
and w ^ T(N, e) with T(7V, e) defined in the formulation of the lemma the function g(x, u) 
monotonically increases in x ^ N, whence it follows that 

sup g(n, u) = lim g(n, u) — lne, ^ u ^ T(N, e), 

n^N n -' 00 

and 

sup r Jt,e--L,n)^ -ee^ 2 / 2 [ " M 2 e^ M = — f e"^ 3 - lV * 2/2 = r 3 (*, e), 
\ \Jn J 2 J 6xV / 
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Q. E. D. 

Finally, the process of computational optimization can be properly organized with the 
help of the following statements. 

LEMMA 5 (Bhattacharya and Ranga Rao, 1976). For any distribution F with zero 
mean and unit variance we have 

p(F, $) ^ sup ( $0) - — X -^] = 0.54093654 . . . 

x>0 V 1 + X J 

LEMMA 6. For any F G JF 3 and n ^ 400 such that (3 3 + 1 ^ 0.1y/n the following 
estimate takes place: 

,(*.,*) < 0.2727- ^ + 2™ 

The statement of lemma 6 is a result of the algorithm described in (Prawitz, 1975b) or 
(Gaponova and Shevtsova, 2009). 

Since the function 

, 1S 0.27276 + 0.2041 , , 

»W = 6TI • 

monotonically increases for k > 0.2041/0.2727 = 0.74. . . and monotonically decreases for 
< k < 0.74, we have 

f 0.2727, fc ^ 0.75, 

supg(fr) = < 

6>i [ 0.4768/(1 + fc < 0.74. 

Thus, from lemma 6 it follows that for all n and /3 3 such that (/3 3 + k)/y/n < 0.05(1 + k) 
inequality (3) holds with C k = 0.2727 for k ^ 0.75 and with C k = 0.4768/(1 + k) for 
k ^ 0.74. In particular, for k = 0.425 we have 

0.3346. 03 + ° r i 426 , if £±°^ < 0.07125. 



The lemmas presented above give the grounds for restricting the domain of the values 
of s — (j3 3 + k)/ y/n by a bounded interval separated from zero (more details will be given 
below) and for looking for the constant C k in the form 

C k = m&xC(e), C{e)=D(e)/e, D(e) = sup {D(e, n) : n ^ n*} , (11) 

where 

D(e,n)= inf D(e - k/y/n,n,t ,T), (12) 

0<i o <l,T>0 

n* = max{l, [(1 + k) 2 /e 2 ]}, 

here is the least integer no less than x. The condition n ^ is a consequence of 
the inequality /3 3 ^ 1. For the estimation of the supremum in n in the definition of D(e), 
lemma 4 is used for N large enough. The computation of the maximum in e is essentially 
based on the property of monotonicity in e of all the functions used for the estimation 
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of \f n {t)\ and r n (t), and hence, on the monotonicity of the function D(e) = eC(e). This 
property makes it possible to estimate max £ C(£:) using the values of C(e) only in a finite 
number of points. In particular, the following statement holds. 

Lemma 7. For all e 2 > £\ > the following inequality is true: 

max C(e) ^ C{e 2 ) ■ — . 



2.2.2 Proof of theorem 1 

Denote 



0.425 p 3 + 0.425 



n \ n 



Then for e ^ 0.07 inequality (6) is a consequence of lemma 6, and for e ^ 1.62 ^ 
0.541/0.335789 it follows from lemma 5. Thus, to compute Ck the maximization with 
respect to e in (11) is conducted on the interval 0.07 ^ e ^ 1.62. To compute the 
supremum with respect to n ^ n„ = [(1.425/e) 2 ] we use lemma 4 with N = 600 for 
e < 0.1, N — 300 for 0.1 < e < 0.2 and N = 100 for e > 0.2. For the mentioned values of 
e we have n*(0.07) = 415, n*(0.1) = 204, n*(0.2) = 51. The maximum with respect to e 
is estimated by lemma 7 and is attained in the two points: n = 5, e ~ 0.822 (/3 3 ~ 1.413, 
t « 0.385, T = 5.755) and n = 8, e w 0.504 (/3 3 = 1, t ~ 0.293, T = 8.911). Both 
extremal values do not exceed 0.335789, whence, theorem 1 is proved. 



2.2.3 Proof of theorem 2 

Denote 



1 (3 3 + l 



n \/n 



Then for e ^ 0.1 inequality (7) is a consequence of lemma 6, and for e ^ 1.78 ^ 
0.541/0.3051 it follows from lemma 5. Thus, to compute Ck the maximization with respect 
to e in (11) is conducted on the interval 0.1 ^ e ^ 1.78. To compute the supremum with 
respect to n ^ n„ = [4/e 2 ] we use the last statement of lemma 4 with N = 200 and 
T(200, e) = min{5.04/e, 3.76/v^}- It turned out, that the extremal value is attained at 
n -> oo and e w 0.985 (t = 0.356, T = 6.147) and it does not exceed 0.3051, Q. E. D. 



3 An improvement of the analog of the Berry-Esseen 
inequality for Poisson random sums 

3.1 The history of the problem 

In this section we will use theorem 1 to improve the analog of the Berry-Esseen inequality 
for Poisson random sums. Let X l5 X 2 ,... be independent identically distributed random 
variables with 

EXi = /x, DX! = (T 2 >0 and E^ 3 = (3 3 < oo. (13) 



12 



Let N x be a random variable with the Poisson distribution with parameter A > 0. Assume 
that for any A > the random variables N\ and Xi, AT 2 , ... are independent. Set 



S\ — X\ + . . . + Xn x 

(for definiteness we assume that S\ = if N\ = 0). Poisson random sums S\ are 
very popular mathematical models of many real objects. In particular, in insurance 
mathematics S\ describes the total claim size under the classical risk process in the 
«dynamical» case. Many examples of applied problems from various fields where Poisson 
random sums are encountered can be found in, say, (Gnedenko and Korolev, 1996) or 
(Bening and Korolev, 2002). 
It is easy to see that 

ES x = \fi, DS X = A(/i 2 + a 2 ). 

The distribution function of the standardized Poisson random sum 

~ = S x - A/i 
A " v/AGu 2 + <r 2 ) 

will be denoted F\(x). 

It is well known that under the conditions on the moments of the random variable X\ 
given above, the so-called Berry-Esseen inequality for Poisson random sums holds: there 
exists an absolute positive constant C < oo such that 

= sup |F>(x) - »(.)| < - t— . (14) 

Inequality (14) has rather an interesting history. Apparently, it was first proved in (Rotar, 
1972a) and was published in (Rotar, 1972b) with C = 2.23 (the dissertation (Rotar, 
1972a) was not published whereas the paper (Rotar, 1972b) does not contain a proof of 
this result). Later, with the use of a traditional technique based on the Esseen smoothing 
inequality this estimate was proved in (von Chossy, Rappl, 1983) with C = 2.21 (the 
authors of this paper declared that C = 3 in the formulation of the corresponding theorem, 
which is, of course, true, but actually in the proof of this theorem they obtained the value 
C = 2.21). 

In the paper (Michel, 1993) the property of infinite divisibility of compound Poisson 
distributions was used to prove that the constant in (14) is the same as that in the 
classical Berry-Esseen inequality. Although Shiganov's estimate C ^ 0.7655 (Shiganov, 
1986), had been known by that time (the original paper by Shiganov had been published 
in Russian even earlier, in 1982), Michel used the previous record value due to Van Beek 
(Van Beek, 1972) and announced in (Michel, 1993) that C ^ 0.8 in (14). Being not aware 
of this paper of Michel, the authors of the paper (Bening, Korolev and Shorgin, 1997) 
used an improved version of the Esseen smoothing inequality and obtained the estimate 
C ^ 1.99. As it has been already noted, the method of the proof used in (Michel, 1993) 
is based on the fact that if for the absolute constant C in the classical Berry-Esseen 
inequality (1) an estimate C ^ M is known, then inequality (14) holds with C = M. 
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This circumstance was also noted by the authors of the paper (Korolev and Shorgin, 1997) 
in which independently of the paper (Michel, 1993) the same result was obtained, but with 
another currently best estimate M = 0.7655. As we noted in Section 1, the best known 
estimate of the absolute constant in the classical Berry-Esseen inequality was obtained 
in (Tyurin, 2009c), (Tyurin, 2009d): C ^ 0.4785. Therefore, following the logics of the 
reasoning used in (Michel, 1993) and (Korolev and Shorgin, 1997) we can conclude that 
inequality (14) holds at least with C = 0.4785. 

In this section we show that actually binding the estimate of the constant C in (14) to 
the estimate of the absolute constant C in the classical Berry-Esseen inequality is more 
loose. Namely, although the best known upper estimate of Co is M = 0.4785 and moreover, 
although the unimprovable lower estimate of Co is ~ 0.4097..., inequality (14) actually 
holds with C = 0.3051. Thus, here we improve the result of (Korolev and Shevtsova, 
2010b) where we proved inequality (14) with C = 0.3450. 

3.2 Auxiliary results 

The following lemma determines the relation between the distributions and moments of 
Poisson random sums and the distributions and moments of sums of a non-random number 
of independent summands. This lemma will be the main tool which we will use to apply 
the results known for the classical case, to Poisson random sums. 

Here and in what follows the symbol = will stand for the coincidence of distributions. 
Also denote v — X/n. 

Lemma 7. The distribution of the Poisson random sum S\ coincides with the 
distribution of the sum of a non-random number n of independent identically distributed 
random variables whatever integer n ^ 1 is: 

X\ + . . . + Xn x = Y v> \ + . . . + Y vn 

where for each n the random variables Y v> i, . . . , Y^ n are independent and identically 
distributed. Moreover, if the random variable X\ satisfies conditions (13), then for the 
moments of the random variable Y v \ the following relations hold: 

E Y v>1 = fiv, D Y V) i = (/i 2 + a 2 )is, 

E\Y u< i -H 3 < +40i/) for n ^ A. 

Proof. The proof is based on the property of infinite divisibility of a compound 
Poisson distribution which implies that for any integer n ^ 1 the characteristic function 
fs x (t) of the Poisson random sum S\ can be represented as 

fo(t) = exp{A(/(f)-l)}= [exp {!/(/(*)-!)}]"= [f Y ^(t)} n , 

where fy vl is the characteristic function of the random variable Y vX . Hence, the 
distribution of each of the summands Y^i, . . . , Y^ n coincides with the distribution of the 
Poisson random sum of the original random variables: 

Y U)k = Xi + ... + X Nv , fc = 1, ... ,n, 
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where N u is the Poisson-distributed random variable with parameter v independent of the 
sequence Xi,X 2 , . . . Hence we directly obtain the relations for the first and the second 
moments of the random variables Y V) \ and X\. Let us prove the relation for the third 
absolute moments. By the formula of total probability we have 

OO 

E \Y Vtl - fiu\ 3 < e- u (i^\n\ 3 + uE \X X - fiu\ 3 + ^ V —E \X X + . . . + X k - H 3 ) • 

k=2 

Consider the second and the third summands on the right-hand side separately. For this 
purpose without loss of generality we will assume that n ^ A, that is, v ^ 1. By virtue of 
the Minkowski inequality we have 



(E \ Xl - H 3 ) 1/3 < (/? 3 ) 1/3 + W= (/? 3 ) 1/3 f 1 + tS 
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Since v ^ 1 and the ratio If^l/iP 3 ) 1 ^ 3 does not exceed 1 by virtue of the Lyapunov 
inequality, we obtain 

E \X X - iiu\ 3 ^ /3 3 (1 + uf ^ /3 3 (1 + 7v). 
To estimate the third summand notice that the Lyapunov inequality yields 

k k 

xA ^ k^ 1 \xi\ r , ij 6 R, % = 1, . . . , k, r ^ 1, 

i=l i=l 

(see, e. g., (Bhattacharya and Ranga Rao, 1976)). With r = 3, this inequality implies 

E \X X + . . . + X k - H 3 < E (|Xi| + ... + \X k \ + |//|i/) 3 ^ 

^ {k + 1) 2 (£;/? 3 + (l^lz/) 3 ) < f3 3 {k + l) 3 
(here we took into account that \fi\ 3 ^ [3 3 and 1/ ^ 1). Thus, 



OO JU 



E \Y Vjl - nv\ 3 ^ z/ 3 |/i| 3 + i/E |Xi - nv\ 3 + — E \X X + . . . + X k - jiu\ 3 < 



fc=2 



^ z^ 3 [l + (8 + i<>] 

where 



k\ 

k=2 



The lemma is proved. 

COROLLARY 1. Under conditions (13) £/ie distribution of the standardized Poisson 
random sum S\ coincides with the distribution of the normalized non-random sum of n 
random variables whatever integer 1 is: 



S\ — —j= y~] z vk 
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, Z v<n are independent and identically 



where for each n the random variables Z v>1 , 
distributed. Moreover, these random variables have zero expectation, unit variance and 
for all n ^ A their third absolute moment satisfies the inequality 

p 3 (l + 40//)^ 



E\Z ul \ d ^ 

' (^2 + ^2)3/2^ 

Proof. According to lemma 7 for any n we have the representation 



(15) 



S\ — A/i 



YvA + • • • + Y v , 



2 ^ a 2 ) 



a 2 )nis 



1 n 
v k=i 



in which the random variables 

Z v ,k 



EY, 



u.k 



D Y, 



u.k 



are independent, identically distributed, have zero expectation and, unit variance. 
Moreover, by virtue of the same lemma for all n ^ A we have the relation 



E|^,i| 



E\Y U +-EY V ^ ^ /3 3 (1 + 40^) /3 3 (l+40^)v^ 



D 1^)3/2 ^ (^ + ^2)3/2^1/2 ^ 2 + ^2)3/2^ 



The corollary is proved. 



3.3 Main result 

THEOREM 2. Under conditions (13) /or any A > we have the inequality 

,„ ^ 0.3051/3 3 

Proof. From lemma 7 and corollary 1 it follows that for any integer n ^ 1 



p(F A , $) = sup 



,fc < £ 



Hence, by theorem 1 for an arbitrary integer n ^ 1 we have 

p(F A ,$K 0.3051^^ + ^. (16) 

Since n ^ 1 is arbitrary, we can assume that n ^ A, making it possible to use estimate 
(15) for the specified n and, in the continuation of (16), to obtain the inequality 

,„ , /3 3 (l + 40A/n) 0.3051 

p(F x , $ ^ 0.3051 — '—^ + — 

^ V ; (/i 2 + a 2 ) 3 / 2 \/A 

Since here n ^ A is arbitrary, letting n->oowe finally obtain 



p(F x , $) < lim 



0.3051 



/? 3 (1 +40A/n) 
(/i 2 + a 2 f/ 2 V\ 



0.3051 



0.3051/3 3 



(/i 2 + a 2 ) 3 / 2 ^' 



Q. E. D. 
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4 Convergence rate estimates in limit theorems for 
mixed compound Poisson distributions 



4.1 Preliminaries 

Let A t be a positive random variable whose distribution depends on some parameter t > 0. 
The distribution function of A t will be denoted G t (x): G t (x) = P(A t < x). By a mixed 
Poisson distribution with a structural distribution G t we will mean the distribution of the 
random variable N(t) which takes values k = 0, 1, ... with probabilities 

oo 

p(jV(t) = fc) = i J e~ x \ k dG t (\), k = 0,1, 2, ... 
o 

Some special examples of mixed Poisson distributions are very well-known. The 
most well-known and most widely used mixed Poisson distribution is, of course, the 
negative binomial distribution (since it was first used in the form of a mixed Poisson 
distribution in (Greenwood and Yule, 1920) to model the frequencies of accidents). 
This distribution is generated by the structural gamma-distribution. Other examples 
of mixed Poisson distributions are the Delaporte distribution with the shifted gamma- 
structural distribution (Delaporte, 1960), the Sichel distribution with the generalized 
inverse Gaussian structural distribution (Holla, 1967), (Sichel, 1971), Willmot, 1987), 
The generalized Waring distribution (Irwin, 1968), (Seal, 1978). The properties of mixed 
Poisson distributions are described in detail in (Grandell, 1997) and (Bening and Korolev, 
2002). 

Let Xi,X 2 ,... be independent identically distributed random variables. Assume that 
the random variables N(t), X 1 , X 2 , ... are independent for each t > 0. Set 

S(t) = X 1 + ... + X m 

(for definiteness we assume that if N(t) = 0, then S(t) = 0). The random variable S(t) 
will be called a mixed Poisson random sum and its distribution will be called compound 
mixed Poisson. 

In what follows we will assume that the random variables Xi,X 2 ,... possess three 
first moments for which we will use the same notation as in Section 3 (see (13)). The 
asymptotic behavior of the distributions of mixed Poisson random sums S(t) when N(t) 
infinitely grows in some sense, is principally different depending on whether /x = or not. 

The convergence in distribution and in probability will be respectively denoted by the 
symbols =>- and 

First consider the case \i = 0. In this case the limit distributions for standardized 
mixed Poisson sums are scale mixtures of normal laws. Without loss of generality, unless 
otherwise indicated, we will assume that a 2 = 1. 

Theorem 3 (Korolev, 1996), (Bening and Korolev, 2002). Assume that A t oo as 
t — > oo. Then, for a positive infinitely increasing function d(t) there exists a distribution 
function H(x) such that 
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if and only if there exists a distribution function G(x) such that for the same function 
d(t) 

G t (xd(t)) G{x) (t->oo) (17) 

and 



H{x) = J $(x/y/y)dG{y), xE 



o 

Now consider the case \i ^ 0. This case is important from the point of view of insurance 
applications. Recall that, in general, DXi = a 2 . Assume that there exist numbers £ G 
(0, oo) and s 6 (0, oo) such that 

EA f = it, DA t = s 2 t, t > 0. (18) 

Then it is easy to make sure that 

ES(t) = fiit, DS{t) = [i{fi 2 + a 2 ) + fi 2 s 2 }t. 

In the book (Bening and Korolev, 2002) a general theorem presenting necessary and 
sufficient conditions for the convergence of compound mixed Poisson distributions was 
proved. The following theorem is a particular case of that result. 

Theorem 4 (Bening and Korolev, 2002). Let ji ^ 0. In addition to the conditions 

p 

on the moments of the structural random variable A t assume that A t — > oo as t — ► oo. 
Then, as t — >• oo 7 compound mixed Poisson distributions converge to the distribution of 
some random variable Z , that is, 

S(t)-fdt ^ 7 
^/[iifx 2 + a 2 ) + fi 2 s 2 }t 

if and only if there exists a random variable V such that 

At -it 



V. 



Furthermore, 




1+ , f S \, n - ; ^ sV ), XER. 

(v 2 + o- 2 )i v^TT^K/ 

It is easy to see that the limit random variable Z admits the representation 



Z = 



1+ » S 



2 C ,2 1-1/2 



(/i 2 + a 2 )t 



x + - » s -v, 

^/{n 2 + a 2 )i + fi 2 s 2 



where X is a random variable with the standard normal distribution independent of V. 

The basic distinctions of the case /i ^ from the case of compound mixed Poisson 
distributions with zero expectations considered above are, first, the necessity of non-trivial 
centering and different normalization required for the existence of non-trivial limit laws 
and, second, the shape of the limit law which in this case has the form of a location 
mixture of normal laws. 
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4.2 Convergence rate estimates in limit theorems for mixed 
compound Poisson distributions with zero mean 

It is easily seen that the distribution of the mixed Poisson random sum S(t) can be 
represented as 



r r Nx \ 

P(S(t) <x)= P(J2 x i < x ) dG tW, x e K. 

n .7—1 



Recall that here we assume that 



EXi = 0, EX 2 = 1, /3 3 = E|Xx| 3 < oo. 



(19) 



(20) 



Let d(t), t > 0, be a positive infinitely increasing function. In this section we will present 
some estimates of the rate of convergence in theorem 3. 
For A > denote 



p(X) = sup 



< X 



Let G(x) be a distribution function such that G(0) = 0. If condition (17) holds, then, 
according to theorem 3, compound mixed Poisson distribution of the mixed Poisson sum 
S(t) normalized by the square root of the function d(t) converges to the scale mixture of 
normal laws in which G(x) is the mixing distribution. Denote 



At = sup 



■(^L< x 



oo 



, St = sup \ G t (d(t)x) — G( 



x) 



THEOREM 5. Assume that conditions (20) hold. Then for any t > we have the 
estimate 

A t < 0.3051 • /3 3 E[A t ]~ 1/2 + 0.5 • 5 t . 

Proof. This statement was first proved in the paper (Gavrilenko and Korolev, 
2006) with a slightly worse constant (also see (Korolev, Bening and Shorgin, 2007). Here 
we present a modified version of the proof. By virtue of representation (19) we have 



A t = sup 



dG(X) 



sup 

x 



sup 

x 



OO , \ OO 

J p(^X s <xy/dmdG t (X)- J S>(^ 

^ ■ ?=1 ' 

tjy \ oo 



_gx 3 <^j dGt (M ( ))j-/*y dG (A) 
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^ sup 

x 



oo 




/ 










Xd(t) 



< 



v/A. 



dG t (\d{t)) 



+ 



sup 



*(^W,(Ad(f))-G(A)] 



Continuing this chain of relations with the use of integration by parts and theorem 2 we 
further obtain 



Aj ^ / sup 



Xd(t) 

7W) Xj 



< x - $fx] 



dG t (\d{t))- 



sup 



[G t (Ad(0) -G(A)]d A $^j 



< J p(X)dG t (\) + sup \G t (Xd(t)) -G(X)\ -sup J d x <&(^=^j 



oo 

^ 0.3051 ■ /3 3 / -^=dG t (A) + 0.5 • sup I G t (Ad(t)) - G(A) I = 
J v A a 



= 0.3051 ■ /5 3 E[A^ 1/2 + 0.5 • <J t , 

Q. E. D. 

As an example of applications of theorem 5 consider the case where for each t > the 
random variable A t has the gamma-distribution. This case is very important in financial 
applications for the asymptotic validation of such popular models of the evolution of 
financial indexes as variance-gamma Levy processes (VG-processes) (Madan and Seneta, 
1990) or two-sided gamma-processes (Carr, Madan and Chang, 1998). 

As is well known, the density of the gamma-distribution with shape parameter r > 
and scale parameter a > has the form 

g r>(T {x) = —— e - ax x r - 1 , x>0. 
1 (r) 

Thus, the mixed Poisson distribution with the mixing gamma-distribution has the 
characteristic function 



e W {y(^ - l^^-e-^-'dy 



a 



T(r) 



exp\ -oi/ l 



- elt ^-y^dy=(u l ~ eU 



a 



a 
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By the re-parametrization 



we finally obtain 



1 — p \ 1 + a J 



which coincides with the characteristic function of the negative binomial distribution with 
parameters r > and p G (0, 1). So, in the case under consideration for each t > the 
random variable N(t) has the negative binomial distribution with parameters r > and 
pe (0,1): 

p(N(t) =n)= C? +n _ lP r (l - p) n , n = 0, 1, 2, . . . . (21) 

Here r > and p G (0, 1) are the parameters and for non-integer r the quantity Cy +n _ 1 is 
defined as 

£m _ T(r + n) 



r+n—l 



n\ ■ T(r) 



In particular, with r = 1, relation (21) determines the geometric distribution. 

The gamma-distribution function with scale parameter a ana shape parameter r will 
be denoted G r>(T (x). It is easy to see that 

G Tt<T {x) = G r>1 {ax). (22) 

The random variable with the distribution function G r>CT (x) will be denoted U(r, a). It 
is well known that 

EU(r,a) = -. 

a 

Fix the parameter r and take U (r, a) as the random variable A t assuming that t = a^ 1 : 

A t = U(r,t~ 1 ). 



As a function d(t) take 
Obviously, we have 



d(t) = EA t = EUir.r 1 ). 



Et/(r,t _1 ) = rt. 
Then with the account of (22) we have 

G t {xd{t)) = P(f/(r,r 1 ) < xrt) = P(U{r, 1) < xr) = P(U(r,r) < x) = G r>r (x). 

Note that the distribution function on the right-hand side of the latter relation does not 
depend on t. Therefore the choice of d(t) specified above trivially guarantees the validity 
of condition (17) of theorem 3. Moreover, in this case 5 t = for all t > 0. 
Now calculate E[Aj] -1 / 2 under the condition 

r > l -. (23) 
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We have 



oo 



-x/t~r-3/2 r( r _ I) 

E[A f ]- 1/2 = E[U(r,t- r )]-^ 2 = I - 1 . dx = 1 2 J 

t r T(r) T(r)y/i 



o 



Thus we obtain the following statement which is actually a particular case of theorem 5. 

Corollary 2. Let the random variable A t have the gamma- distribution with shape 
parameter r > and scale parameter a = 1/t, t > 0. Assume that conditions (20) and 
(23) hold. Then for each t > we have 



sup 



oo 

P(S(t) < xVri) - J $ (^-J dG 

^ 



•(!/) 



< 0.3051 ^— 21 ■ -p. 

r(r) Vt 



If r = 1, then the random variable 

N(t) = N 1 (U(l,r 1 )), t>0, 

has the geometric distribution with parameter p = As this is so, the limit (as t — > oo) 
distribution function of the standardized geometric sum S(t) is the Laplace distribution 
with the density 

l(x) = -L e -^, x e R. 

The distribution function corresponding to the density l(x) will be denoted L(x), 

\e^ x , if x < 0, 

1 - ie-^, if x > 0. 
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Corollary 3. Let £/te random variable A t have the exponential distribution with 
parameter a = 1/t, t > 0. Assume that conditions (20) /io/d Then for each t > 

sup |P(S(t) < xv^) - L(x)\ < 0.5408 ■ 

5 Convergence rate estimates in limit theorems for 
mixed compound Poisson distributions with non- 
zero mean 

Here we will present some estimates of the rate of convergence in theorem 4. 

5.1 The case of structural random variables with finite variance 

Under assumptions (18) denote 

F t (x) = P\ - <x 

y/[£{fM 2 + a 2 ) + fJ, 2 S 2 ]t 
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Pt = sup 



,2o2 



F t (x)-E<S>lxJl + — ^ . 



G*(u) = P(y < u), 5 t = sup G t (ysVi + et) -G*( 



Theorem 6. Lei EXi = // ^ 0, DXi = cr 2 , E|Xi| 3 = (3 3 < oo, E |V| < oo. T/ien for 
any t > we have 



Pt < S t + -p • inf 



0.3051/3 3 



where 



Q(e) = max < - 



e' (1 + vT^i) ^2^(1 - e) J ' 



Proof. A similar statement with slightly worse constants was first proved in the 
paper (Artyukhov and Korolev, 2008). Here we present a modified version of the proof. As 
above, let N\ be a random variable with the Poisson distribution with parameter A > 
independent of the sequence X\, X2, ■ ■ ■ Then we can write 



Pt = sup 



S{t) - fdt 



sup 



v/[% 2 + a 2 ) + fi 2 s 2 ]t 



< x } - E$ \x\ I 



\i 2 s 2 



fisV 



(/i 2 + cr 2 )£ v/(/i 2 + a 2 )£ 



^[£{fi 2 + a 2 )+fi 2 s 2 }t 



< x dG t (X)- 



E$ x* I 



H 2 s 2 



fisV 



Fix an arbitrary e G (0, 1). Then 



Pt = sup 

x 



Sn x — p£t 



\<(i-e)et 



^/[i(fi 2 + a 2 ) + fi 2 s 2 }t 



< x dGt(X)- 



+ 



Sn x — p-£t 



\>{l+e)U 



^[£(fi 2 + a 2 ) + fi 2 s 2 ]t 



Sn x — i^£t 



(l-e)£t^A^(l+e)a 



^[£(fi 2 + a 2 ) + fi 2 s 2 }t 



< x dG t {X)+ 



< x dG t (X)- 



E$ xjl 



H 2 s 2 



[isV 



(/i 2 + o- 2 )£ y/QF+jSfi 
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^ sup 



S Nx — jj.it 



\<{l-e)£t 



^[i(fi 2 + a 2 ) + 

Sn x — V-it 



< x dG t (X)+ 



A>(l+e)tt 



+ sup 

x 



< x dG t (X) 



+ 



(l-e)ft^A^(l+e)£f 



^/[i(fi 2 + a 2 )+fi 2 s 2 }t 



-E$( 



374/1 



H 2 s 2 



< x dGt(X)- 



[isV 



V V ^ 2 + ° 2 Y V^ 2 + ° 2 Y 



< p 



it 



>e + 



sup 



S Nx - fiit 



(l-e)ft^A^(l+e)ft 



^(/Z 2 + <7 2 ) + //V]t 



< x dG t {X)- 



E$ x, 1 



H 2 s 2 



fisV 



Further, 



sup 



G" 2 + ° 2 Y v/(/i 2 + a 2 )i 
Sn x — f-it 



(24) 



(l-e)«^A^(l+e)ft 



^[i^ 2 + a 2 )+fi 2 s 2 ]t 
-E$( 



V V (/^ 2 + ^ 2 K vV + ^ 2 K 



^ sup 

a; 



(l-e)«<A<(l+e)ft 



"V 1 

^ N\ — ^X 



/i 2 S 2 



< x ) dG t (X)- 
fisV 



Vxjjj, 



2 ^a 2 ) 



< 



A/y I 



H 2 s 2 



[IS 



X-it 



+ sup 

x 



i/yl 



\i 2 s 2 



(v 2 + ° 2 y v^ 2 + ° 2 y s ^ 

fis X — it 



A 



/i 2 S 2 



dG t (A)+ 
X-it 



(l-e)ft^A^(l+e)ft 



E$( xWl 



/i 2 S 2 



(/i 2 + <J 2 K + a 2 )£ sV* 



dG t (A)- 



V V (/^ 2 + ^ 2 K v/(/i 2 + a 2 )£ 



/i + / 2 
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Consider I\. Denote 



y?s 2 /xs X — it 

X\ X \I + (/i 2 + a 2 )i ~ 7p + a 2 )i ' ~s^t 



Then ii can be rewritten in the form 



l\ = sup 

y 



Sn x — fiX 

v / a(^+^ 2 : 



< y - $(y) 



sup 

y 



dG t (X) ^ 



dG t (X). 



(l-e)&<A<(l+e)& 

To estimate the integrand on the right-hand side of the latter inequality we use theorem 
2 and obtain 



0.3051/3 3 
(/i 2 + a 2 ) 3 / 2 



0.305L3 3 



X£(l-e)tt 



(n 2 + a 2 ) 3 / 2 ^/(l-e)£t 



(25) 



Consider J 2 . We have 



I 2 ^ sup 



(l-e)£t^A^(l+e)ft 



\i 2 s 2 



(IS 



X-tt 



(/i 2 + vV + a 2 )i sy/i 



H 2 s 2 



[IS 



x-et 



+ sup 



(/i 2 + ^ 2 K vV + ° 2 Y s ^ 



$ \XA 1 



jJ 2 S 2 



dG t (X)- 
X-lt 



(l-e)£t^A^(l+e)ft 



/i 2 S 2 



(/u 2 + a 2 K v^T^K 



dG t (A)- 



(/i 2 + ^ 2 )£ v^TT^ 



/21 + -^22- 



Denote 



Then 



z = X\ 1 



/i 2 s 2 



{IS 



X-U 



(/i 2 + vV + ^ ' 



hi ^ 



sup 



dG t (A). 



(26) 



(l-e)ft<A^(l+e)ft 

Consider the integrand in (26) . By the Lagrange formula we have 



A" 1 



tp[ez + (i-9)z\lj 



(27) 
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for some 6 G [0, 1] where 



ip(x) = &(x) 



-x 2 /1 



is the standard normal density. The function <f(x) = <fi(\x\) monotonically decreases as 
|x| increases. Therefore the function ip on the right-hand side of (27) attains its maximum 
value in 9 G [0, 1] at that value of its argument, whose absolute value is minimum. But 
the argument of the function (p on the right-hand side of (27) is itself a linear function of 
9. Therefore, the minimum absolute value of this argument is attained either at 9 = or 
at 9 = 1. But at 9 = 1 we have 



9 1 



9 1 



while at 9 = we have 



In the definition of I21 A satisfies the inequality A ^ (1 + e)£t. Therefore, 

[E 1 

X ' y/T+e' 



Hence, 



— ^ min < 1 



VT 



VT 



Therefore in J 2 i we have (see (27)) 



sup 



Furthermore, 



vrn 



■ sup \z\ip 



1 - 



VT 



Therefore the supremum in (28) is attained at z = ±\/T+~e and equals 



sup \z\ip 



VTT 



2vre 



Thus, 



hi < 



2vre 



(l-e)«^A^(l+e)« 



2ire 



(l-e)«^A^(l+e)fl 



Va 



Va 



dG t (A) 

Ji + Vx 



£t + VX 



dG t (X) 



(28) 
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i + e r \x-u\ 

27re J ^/xl^lt + y/X) 

(l-e)ft<A<(l+e)ft ' 



dG t (X) < 



1 + e 



27re V /T^7(l + V /T^) 





A 


/ 


S" 1 



1 + e 



27re(l - e) 1 + ^l^~e 

Here the symbol 1(A) denotes the indicator of a set A. 
Consider J 22 . We have 



)&<A<(1- 






At 

— - 1 

ft 




A t 

— - 1 

ft 



dG t (X) 



< e 



(2 



\i 2 s 2 



[IS 



A:|A-l|<e 



dG t (A) 



[i 2 s 2 



[IS 



X-lt 



dG t (X) 



X:\^KeiVt/s 



dG t (vsVi + £t). 



[i 2 s 2 



pLSV 



Therefore, 



I 22 = SUp 



$ hrWl 



H 2 s 2 



jisv 



\v\^.tt\/l/s 



(^ 2 + ° 2 Y VW+^Yi 



dG t (vsVi + £t)- 



$ \x* 1 



H 2 s 2 



jJSV 



dG*(v) < 



sup 



[dG t (vsVi + it) -G* 



H 2 s 2 



fisv 



\v\^e£y/i/s 



(/i 2 + ° 2 Y V^ 2 + a 2 )i 



+ 



+ sup 



$ \x* 1 



\i 2 s 2 



fISV 



\v\>eiVi/s 

By integration by parts we obtain 



(/i 2 + * 2 )* + ° 2 Y 



dG*(v) 



-^221 + -^222- 



J 22 i ^ sup 



J [G t (vsVi + £t) -G*(v)]d v ® IxJl 



\i 2 s 2 



/jLSV 



\v\^el^t/s 



(/i 2 + ° 2 Y vV + a 2 )i 
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^ sup I G t (vsVt + it) - G*(v)\ = St. (30) 

V 

Note that in (24) we can apply the Markov inequality and obtain the estimate 











i.E 


i — 1 


e 









> e 



(31) 



Further, again applying the Markov inequality we can make sure that 

sE\V\ 



/ 222 < P(\V\>ef.Vi/s) <: 



(32) 



Now unifying (24), (25), (26), (29), (30), (31) and (32) we finally obtain that for any 
e G (0, 1) there holds the inequality 



pt < St 



0.3051/3 3 



-If- 

yft\{^ + a 2 fl 2 yj{l-e)i d 



■ BV\ 



+ E 



A t 

— - 1 

et 



max < - 



1 



1 + e 



e'V27re(l-e) 1 + JT^l J ' 



whence we obviously obtain the statement of the theorem since the Lyapunov inequality 
obviously implies that for each t > 



A* 

et 



iyft 



At -it 



eVt 



y/t J £y/t 



The theorem is proved. 

If we additionally assume that the family of random variables 

At -it 



iy/t 



t>o 



is uniformly integrable, then by the Lyapunov inequality we obtain the inequality 

At -it 



E\V\ = lim E 

t— *oo 



< Um WD 1^-^- 

W V Sy/t 



A 



1. 



(33) 



Hence, from theorem 6 we obtain the following result. 

COROLLARY 4. In addition to the conditions of theorem 6, let (33) hold. Then for any 
t > 

1 , ( 0.305L3 3 8/1 



Pt < s t 



inf 



y/i ,6(0,1) I (/i 2 + a 2 ) 3 /V(l 1 



Q« 



where 



Q(e) = max < - 



y/T+e 



e' (1 + vT^i) y/2ne(l-e) 
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5.2 The case of structural random variables with infinite variance 



Assumption (18) which guarantees the existence of the variance of the structural random 
variable At is not crucial. An analog of theorem 6 can be proved for the case where only 
the existence of the mathematical expectation of A t is assumed. Namely, the following 
theorem holds. 

p 

THEOREM 7. Let /i ^ 0. Assume that EA t = t and A t — > oo as t — »• oo. Then, 
as t — > oo, the distributions of normalized mixed Poisson random sums converge to the 
distribution of some random variable Z , that is, 

S(t) - u „ 

if and only if there exists a random variable V such that 

A t -t 



Vi 



v. 



Moreover, 

P(Z < x) = E$ 
Relation (34) means that in theorem 8 



x — pV 



x G 



(34) 



Z = yZ/i 2 + a 2 -X + fiV 



where the random variables X and V are independent and X has the standard normal 
distribution. 



By analogy with the notation introduced above, denote 



Pt = sup 



S(t) - pt 



< X 



x — pV 



a 2 + ji 2 



8 t = sup 



G t {vVi + t))-G*(v) 



Theorem 8. Assume that /3 3 < oo, EA t = t,t>0, and E \V\ < oo. Then 



Pt < St 



inf 



0.3051/3 3 



E|y| 



y/t 66(0,1) i (p 2 + a 2 )V 2 VT^l 



Q(e)E 



At-t 



Vi 



where 



Q(e) = max < - 



Vi + 



e (1 + VT^e) V 2 ^ 1 ~ e ) 
The proof of theorem 8 differs from the proof of theorem 6 only in notation. 
As an example of the situation in which theorems 7 and 8 are valid, but theorems 4 



and 6 are not, consider the case where 

At = max{0, VtV + t} 



1 /2a + 1 



2t a / 2 



a 



Vt-l 
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with 2 < a < 3 and V being the random variable with the density 



yK 1 2(\x\ + l) a 



It can be easily verified that EA t = t for any t > 0, but the second moment of A t is 
infinite due to that the second moment of the random variable V does not exist (and 
hence, the second moment of the mixed Poisson random sum S(t) with the structural 
random variable A t does not exist) . However, it can be easily seen that 



as t — > oo. This case is an illustrative example of an interesting and non-trivial fact: 
unlike the classical summation theory, for sums with a random number of summands (in 
particular, for mixed Poisson random sums) with infinite variances the existence of non- 
trivial weak limits is possible under the normalization of order t 1 / 2 which is «standard» 
in the classical theory only for sums with finite variances. 

The authors have the pleasure to express their gratitude to Margarita Gaponova who 
carried out the supplementary computations resulting in lemma 6. 
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